home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dorghr.z / dorghr
Encoding:
Text File  |  2002-10-03  |  4.0 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDOOOORRRRGGGGHHHHRRRR((((3333SSSS))))                                                          DDDDOOOORRRRGGGGHHHHRRRR((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DORGHR - generate a real orthogonal matrix Q which is defined as the
  10.      product of IHI-ILO elementary reflectors of order N, as returned by
  11.      DGEHRD
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DORGHR( N, ILO, IHI, A, LDA, TAU, WORK, LWORK, INFO )
  15.  
  16.          INTEGER        IHI, ILO, INFO, LDA, LWORK, N
  17.  
  18.          DOUBLE         PRECISION A( LDA, * ), TAU( * ), WORK( * )
  19.  
  20. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  21.      These routines are part of the SCSL Scientific Library and can be loaded
  22.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  23.      directs the linker to use the multi-processor version of the library.
  24.  
  25.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  26.      4 bytes (32 bits). Another version of SCSL is available in which integers
  27.      are 8 bytes (64 bits).  This version allows the user access to larger
  28.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  29.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  30.      only one of the two versions; 4-byte integer and 8-byte integer library
  31.      calls cannot be mixed.
  32.  
  33. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  34.      DORGHR generates a real orthogonal matrix Q which is defined as the
  35.      product of IHI-ILO elementary reflectors of order N, as returned by
  36.      DGEHRD: Q = H(ilo) H(ilo+1) . . . H(ihi-1).
  37.  
  38.  
  39. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  40.      N       (input) INTEGER
  41.              The order of the matrix Q. N >= 0.
  42.  
  43.      ILO     (input) INTEGER
  44.              IHI     (input) INTEGER ILO and IHI must have the same values as
  45.              in the previous call of DGEHRD. Q is equal to the unit matrix
  46.              except in the submatrix Q(ilo+1:ihi,ilo+1:ihi).  1 <= ILO <= IHI
  47.              <= N, if N > 0; ILO=1 and IHI=0, if N=0.
  48.  
  49.      A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
  50.              On entry, the vectors which define the elementary reflectors, as
  51.              returned by DGEHRD.  On exit, the N-by-N orthogonal matrix Q.
  52.  
  53.      LDA     (input) INTEGER
  54.              The leading dimension of the array A. LDA >= max(1,N).
  55.  
  56.      TAU     (input) DOUBLE PRECISION array, dimension (N-1)
  57.              TAU(i) must contain the scalar factor of the elementary reflector
  58.              H(i), as returned by DGEHRD.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDOOOORRRRGGGGHHHHRRRR((((3333SSSS))))                                                          DDDDOOOORRRRGGGGHHHHRRRR((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      WORK    (workspace/output) DOUBLE PRECISION array, dimension (LWORK)
  75.              On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
  76.  
  77.      LWORK   (input) INTEGER
  78.              The dimension of the array WORK. LWORK >= IHI-ILO.  For optimum
  79.              performance LWORK >= (IHI-ILO)*NB, where NB is the optimal
  80.              blocksize.
  81.  
  82.              If LWORK = -1, then a workspace query is assumed; the routine
  83.              only calculates the optimal size of the WORK array, returns this
  84.              value as the first entry of the WORK array, and no error message
  85.              related to LWORK is issued by XERBLA.
  86.  
  87.      INFO    (output) INTEGER
  88.              = 0:  successful exit
  89.              < 0:  if INFO = -i, the i-th argument had an illegal value
  90.  
  91. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  92.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  93.  
  94.      This man page is available only online.
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.